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reported.  Damage  Initiation  site  and  mode  are  compared  with  the  analytically 
predicted  three  dimensional  stress  field  at  the  edge  of  the  hole.  The  location 
of  first  observable  damage  in  all  laminates  was  an  experimentally  repeatable 
layup  characteristic.  All  laminate  types  showed  a  nearly  horizontal  initial 
region  of  damage  growth  followed  by  a  second  region  of  gradually  increasing 
damage  area  growth  with  loading  cycles,  and  a  final  region  of  very  rapid  growth 
just  prior  to  failure.  Post  fatigue  inspection  showed  damage  concentrated  mainl 
in  the  angle  plys.  No  well  defined  relationship  could  be  established  however, 
between  specimen  fracture  and  final  damage  area.  Life  to  first  damage  and  life 
to  specimen  failure  for  the  four  laminate  types  showed  that  both  stacking 
sequence  and  ply  orientation  had  an  affect  on  the  response  of  a  specimen  to 
cyclic  compression  loads. 


SUMMARY  OF  MAJOR  FINDINGS 


Fracture  of  [02/±6]g  Laminates  with  a  6.35  mm  Diameter  Hole 


1 .1  Eight  ply  [O2/±0]g  laminates  were  fabricated  into  tensile  coupons 
each  containing  a  6.35  mm  diameter  hole.  The  cross-ply  angles, 

6,  were  15,  30,  45,  60,  75,  and  90  degrees.  There  were  five 
specimens  for  each  angle  6.  Figure  1  shows  the  nominal  dimen¬ 
sions  of  the  coupon  type  of  specimen  used  in  these  investigations. 

1 .2  The  ultimate  stresses  are  shown  in  Figure  2  and  the  ultimate 
strains  measured  at  fracture  are  shown  in  Figure  3. 

1 . 3  It  is  seen  from  Figures  2  and  3  that  the  laminates  [02/^30]^ 
exhibit  two  distinct  levels  of  behavior.  One  group  (3  specimens) 
failed  at  an  appreciably  lower  fracture  stress  and  ultimate  strain 
than  the  other  group  (2  specimens).  As  can  be  seen,  the  ultimate 
stresses  versus  6  are  aligned  into  two  distinct  curves  with  a 
marked  discontinuity  at  6=30®.  The  ultimate  strains  versus  6 

are  similarly  aligned. 


1 .4  The  stress-strain  curves  for  the  specimens  with  6  less  than  30° 
are  essentially  linear  to  failure  whereas  for  the  specimens  with 
0  greater  than  30°  the  stress-strain  curves  are  essentially 
bilinear  (see  Figure  4)  to  failure.  For  the  specimens  of 
laminate  [02/±30]^ ,  the  three  specimens  which  failed  at  higher 
stress  exhibited  the  bilinear  form  of  stress-strain  behavior 


while  the  other  two,  which  failed  at  lower  stress,  were 


linear  to  failure. 
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Crack  Growth  of  Splits  in  Uni -directional  Laminates  with  a  Hole 

2.1  In  an  earlier  investigation  the  splitting  of  uni-directional 
laminates  containing  a  hole  caused  by  the  static  (one-time) 
application  of  tensile  loads  was  investigated.  In  this  research 
4-ply  uni-directional  laminates  containing  a  3.2  mm  (.125  inch) 
diameter  hole  were  cyclically  loaded  at  constant  amplitude  to 
determine  the  growth  of  the  split  with  the  number  of  cycles. 
Seventeen  specimens  were  cycled  between  8%  and  80%  of  the  nominal 
splitting  stress  of  330  MPa.  Three  specimens  were  cycled  between 
7%  and  70%,  and  four  were  cycled  between  6%  and  60%.  Most  tests 
were  terminated  at  50,000  cycles  with  a  few  being  continued  to 
500,000  cycles. 

2.2  Typical  crack  length  vs.  cycle  data  is  shown  in  Figure  5.  The 
splits  do  not  initiate  symmetrically  from  the  edge  of  the  hole 
and  as  indicated  in  Figure  5,  the  splits  are  located  relative 
to  the  hole  in  four  locations.  The  curve  labeled  T  is  the  sum 
of  the  lengths  labeled  1  and  2  where  B  is  the  sum  of  3  and  4. 

2 ■ 3  The  crack  growth  data  has  been  fitted  to  the  equation 

2a  =  A  log  N-K 

where 

2a  is  the  total  crack  length 
N  is  the  number  of  cycles 


'A'  and  ' K'  are  constants  that  depend  upon  the  cyclic  stress 
levels.  This  equation  was  used  to  extrapolate  backwards  to  Ng, 
the  number  of  cycles  for  crack  initiation. 

2.4  The  values  of  F,  K  and  Ng  for  the  three  cyclic  load  levels  are 
shown  in  Table  IV.  It  should  be  remembered  that  the  number  of 
specimens  for  two  smallest  load  levels  are  quite  small. 

III.  Static  Tensile  Behavior  of  [0/16]^  and  [±6/0]^  Laminates  with  Holes 

3.1  This  research  experimentally  determined  the  fracture  behavior 
of  specimens  containing  various  sized  holes  in  the  laminate 
systems  [0/+e]j  and  [±9/0]^  with  6=0,  5,  10  and  15  degrees. 
Tables  V  and  VI  summarize  the  testing  program. 

3.2  The  fracture,  i.e.  failure,  stresses  for  all  the  specimens  are 
summarized  in  Tables  VII  through  XVI. 

3 . 3  The  effect  of  hole  size  on  the  fracture  strength  has  been 
correlated  against  the  relation 

0^  =  H|.(2a)''" 

where 

2a  is  the  hole  diameter  in  millimeters 

is  a  "composite  toughness"  with  dimensions  MPa(mm)^ 

and  m  is  a  parameter  analogous  to  the  factor  of  .5  (the 
singularity)  in  linear  elastic  fracture  mechanics. 
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Table  XVI  summarizes  the  value  of  and  m  obtained  in  these 
investigations. 

IV.  Damage  Initiation  and  Propagation  During  Compressive  Fatigue  of  Flawed 

Graphite/Epoxy  Composites 

4.1  Compression-compression  fatigue  tests  of  four  different  graphite/ 
epoxy  layups,  (±45/0)^,  (0/±45)^,  (±30/0)^  and  (0/±30)^,  have 
been  conducted.  Buckling  of  these  axially  loaded  laminates  was 
prevented  by  the  use  of  a  sandwich  specimen  composed  of  a  2.5  cm. 
thick  aluminum  honeycomb  core  faced  with  two  5x30.5  cm.  G/E 
coupons.  Damage  initiation  and  propagation  near  the  1.25  cm. 
diameter  hole  in  each  G/E  face  sheet  was  monitored  and  recorded 
during  testing.  Inspection  of  the  laminates  was  by  Moire  inter¬ 
ferometry,  a  method  of  optically  producing  a  visible  contour 

map  of  the  small  out-of-plane  surface  displacements  resulting 
from  subsurface  delamination.  Specimen  cross  sections  were 
microscopically  examined  to  compare  the  observed  Moire  patterns 
with  the  actual  damage  developed  during  cyclic  loading. 

4.2  A  comprehensive  set  of  damage  progression  photographs  are 
presented.  These  pictures,  along  with  plots  of  damage  area 
growth,  show  that  each  of  the  four  laminate  types  displays  its 
own  characteristic  pattern  of  initiation  and  propagation.  No 
well  defined  relationship  between  specimen  fracture  and  final 
damage  area  was  observed,  however,  the  rapid  growth  of  damage 
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FIG.  1  ;  THE  PHYSICAL  CHARACTERISTICS  OF  THE 
COUPON  SPECIMEN 
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TABLE  I 

TENSILE  FRACTURE  SPECIMENS 


LAMINATE 

HOLE  DIAMETERS 
mm 

TOTAL  NO. 

OF  SPECIMENS 

-- 

-- 

[0]6 

-- 

10 

[034 

-- 

5 

[02/±9]5 

6=0,15,30,45,75,90 

6.35 

33 

[o/±s]g 

6=0,5 , 10,15 

0,3.18,6.35,9.53,12.70 

190 

[^e/0]^ 

6=0,5,10,15 

0,3.18,6.35,9.53,12.70 

75 
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TABLE  II 

TENSION-TENSION  CYCLICALLY  LOADED  SPECIMENS 
[0]^  with  3.2  mm  holes 


LOADING  NO.  OF  SPECIMENS 

Q%  to  80%  of  splitting  stress  17 

7%  to  70%  of  splitting  stress  3 

6%  to  60%  of  splitting  stress  4 


splitting  stress  =  333  MPa 
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TABLE  in 

COMPRESSION-COMPRESSION  CYCLICALLY  LOADED  SPECIMENS 
Specimens  with  12.7  mm  Diam.  Hole 


LAMINATE 

NO.  OF  SPECIMENS* 

[±45/0]^ 

CYCLED  TO  FAILURE 

22 

[±45/0]^ 

CYCLED/RESIDUAL  STRENGTH 

26 

C0/±45]^ 

CYCLED  TO  FAILURE 

3 

[0/±45]^ 

CYCLED/SECTIONED 

1 

[±30/0]g 

CYCLED/SECTIONED 

2 

[0/130]^ 

CYCLED/SECTION 

2 

[0/±30]^ 

CYCLED/RESIDUAL  STRENGTH 

1 

*Each  specimen  has  two  laminates. 
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TABLE  IV 

PARAMETERS  FOR  GROWTH  OF  SPLITS  UNDER  CYCLIC  LOADS 


CYCLIC  LOAD  LEVEL 

NO.  OF  SPECIMENS 

4 

K 

NO 

8%  to  80% 

17 

15 

28 

65 

1%  to  70% 

3 

9 

26 

600 

6%  to  60% 

4 

5 

20 

3000 

TABLE  V 


TESTING  PROGRAi^  USED  FOR  IIP/SSTIGATION 


Specimen 

N.  Type 

Nominal 

Hole  Size 

1  '  • 

(0)5 

■ 

(0/i5)3 

(o/iio)^ 

(0/il5)^ 

w 

No  Hole 

U6-X 

5B1-X 

1031-X 

1531-X 

3.175  irun 

U6-H2-X 

531-H2-X 

1031-H2-X 

15B1-K2-X 

6.35  mm 

u6-:-:4-x 

1031-H4-X 

9.525  mm 

U6-H6-X 

5B1-H6-X 

1031-H6-J. 

1531-H6-X 

12.79  mm 

, 

U6-HS-X 

1031-EB-X 

15B1-K8-X 

X  =  Specimen  numcer,  1-10 

10  specimens  -were  -ade  for  each  group  (each  "box) 
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TABLE  VI 


X  =  Speciner.  number,  I-5 

5  specimens  were  made  for  each  group  (each  bon) 


,E  STRESSES  OP  UNPLAV/ED  (0)^,  (0/15 (0/+10)^,  AND  (0/+15)„  SPECIMENS 
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TABLE  VIII 


FAILURE  STRESSES  0?  UUF 
(±5/0)g,  (±io/o)g,  (115/0)3 


LAV.LD 

SPECIMEN 


Specimen 

l-U  (^:Pa) 

5AI-I 

1104 

5A1-2 

1232 

5A1-3 

OQO 

^  y 

5A1-4 

12SS 

5A1-5 

1048 

Average 

1  1134  t  122 

ICAl  1 

10^6 

lCAl-2 

cpc; 

lQAl-3 

5^1 

lOAl-4 

lOAl- 

i 

1103 

Average 

1C16  +  57^ 

15A1-1 

541 

15A1-2 

871 

15A1-3 

578 

15A1-4 

737' 

15A1-5 

830 

Average 

851  7"^ 
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TABLE  IX 

FAILURE  STRESSES  OF  (0)^  SPECIMENS  WITH  HOLES 


Specimen 

o^(MPa ) 

Specimen 

o^(MPa) 

U6-H2-1 

1661 .8 

U6-H6-1 

1546.9 

U6-H2-2 

1576.5 

U6-H6-2 

1752.6 

U6-H2-3 

1  605.6 

U6-H6-3 

1  354.3 

U6-H2-4 

1650.7 

U6-H6-4 

1553.6 

U6-H2-5 

1670.7 

U6-H6-5 

1501  .9 

U6-H2-6 

1635.7 

U6-H6-6 

1366.6 

U6-H2-7 

1551  .8 

U6-H6-7 

1585.8 

U6-H2-8 

1699.6 

U6-H6-8 

1719.3 

U6-H2-9 

1440.3 

U6-H6-9 

1561  .9 

U6-H2-10 

1652.1 

U6-H6-10 

1637.9 

Average 

1614.5 

Average 

1558.1 

Std,  Dev. 

±75.7 

Std.  Dev. 

±130.1 

U6-H4-1 

1 777.3 

U6-H8-1 

1352.7 

U6-H4-2 

1446.7 

U6-H8-2 

1640.0 

U6-H4-3 

1  465.6 

U6-H8-3 

1461  .7 

U6-H4-4 

1446.9 

U6-H8-4 

1445.3 

U6-H4-5 

1  395.6 

U6-H7-5 

1464.0 

U6-H4-6 

1313.0 

U6-H8-6 

1618.0 

U6-H4.7 

1723.9 

U6-H8-7 

1613.4 

U6-H4-8 

1561  .9 

U6-H8-8 

1644.9 

U6-H4-9 

1483.3 

U6-H8-9 

1538.0 

U6-H4-10 

1513.1 

U6-H8-10 

1257.4 

Average 

1512.7 

Average 

1503.5 

Std.  Dev. 

±142.4 

Std.  Dev. 

±131  .2 
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TABLE  X 

FAILURE  STRESSES  OF  (0/i5)^  SPECIMENS  WITH  HOLES 


Specimen 

a^(MPa) 

Specimen 

a.p(MPA) 

5B1-H2-1 

1014.7 

5B1-H6-1 

762.4 

5B1-H2-2 

987.6 

5B1-H6-2 

725.6 

5B1-H2-3 

984.1 

5B1-H6-3 

747.9 

5B1-H2-4 

1015.9 

5B1-H6-4 

729.1 

5B1-H2-5 

1020.2 

5B1-H6-5 

727.2 

5BI-H2-6 

922.1 

5B1-H6-6 

805.3 

5B1-H2-7 

984.1 

5B1-H6-7 

698.9 

5B1-H2-8 

986.0 

5B1-H6-8 

763.0 

5B1-H2-9 

1049.0 

5B1-H6-9 

722.5 

5B1-H2-10 

921  .1 

5B1-H6-10 

757.9 

Average 

Std.  Dev. 

988.5 

±40.9 

Average 

Std.  Dev. 

743.9 

±29.8 

5B1-H4-1 

91  7.7 

5B1-H8-1 

666.1 

5B1-H4-2 

811  .1 

5B1-H8-2 

695.5 

5B1-H4-3 

811  .5 

5B1-H8-3 

630.6 

5B1-H4-4 

824.7 

5B1-H8-4 

704.7 

5BI-H4-5 

1020.5 

5B1-H8-5 

712.4 

5B1-H4-6 

810.0 

5B1-H8-6 

719.0 

5BI-H4-7 

850.9 

5B1-H8-7 

639.9 

1 

5B1-H4-8 

826.2 

5B1-H8-8 

666.8 

5BI-H4-9 

905.3 

5B1-H8-9 

735.4 

5BI-H4-10 

840.6 

5B1-H8-10 

700.2 

Average 

Std.  Dev. 

861  .9 
±67.7 

Average 

Std.  Dev. 

687.1 

±34.7 
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TABLE  XI 

FAILURE  STRESSES  OF  (0/±10)^  SPECIMENS  WITH  HOLES 


Specimen 

o^(MPa ) 

Specimen 

a^(MPa ) 

10B1-H2-1 

696.6 

10B1-H6-1 

581  .7 

10B1-H2-2 

828.7 

10B1-H6-2 

617.7 

10B1-H2-3 

831  .2 

10B1-H6-3 

543.2 

10B1-H2-4 

707.6 

10B1-H6-4 

635.9 

10B1-H2-5 

722.6 

10B1-H6-5 

580.0 

10B1-H2-6 

762.3 

10B1-H6-6 

595.3 

10BI-H2-7 

753.8 

10B1-H6-7 

526.9 

10B1-H2-8 

787.9 

10B1-H6-8 

559.6 

10B1-H2-9 

885.2 

10B1-H6-9 

630.8 

1031-H2-10 

763.9 

10B1-H6-10 

511  .3 

Average 

774.0 

Average 

578.2 

Std.  Dev. 

±60.1 

Std.  Dev. 

±43.0 

10B1-H4-1 

739.5 

10B1-H8-1 

519.4 

10B1-H4-2 

650.7 

10B1-H8-2 

574.4 

10B1-H4-3 

703.8 

10B1-H8-3 

493.7 

10B1-H4-4 

669.5 

10B1-H8-4 

497.9 

10BI-H4-5 

707.7 

10B1-H8-5 

592.6 

10B1-H4-6 

670.8 

10B1-H8-6 

579.8 

10B1-H4-7 

607.9 

10B1-H8-7 

574.1 

10B1-H4-8 

803.8 

10B1-H8-8 

466.3 

10B1-H4-9 

635.3 

10B1-H8-9 

510.7 

10B1-H4-10 

697.6 

10B1-H8-10 

593.3 

Average 

688.7 

Average 

540.2 

Std.  Dev. 

±55.8 

Std.  Dev. 

±47.3 
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TABLE  XII 

FAILURE  STRESSES  OF  (0/:15)^  SPECIMENS  WITH  HOLES 


Specimen 

a^(MPa) 

Specimen 

c^(MPa) 

15B1-H2-1 

613.00 

15B1-H6-1 

414.70 

15B1-H2-2 

600.20 

15B1-H6-2 

453.80 

15B1-H2-3 

681 .90 

15B1-H6-3 

474.60 

15B1-H2-4 

544.60 

15B1-H6-4 

404.20 

15B1-H2-5 

635.10 

15B1-H6-5 

440.30 

15B1-H2-6 

650.00 

15B1-H6-6 

467.80 

15B1-H2-7 

670.80 

15B1-H6-7 

489.80 

15BI-H2-8 

711 .50 

15B1-H6-8 

451 .90 

15B1-H2-9 

702.30 

15B1-H6-9 

391  .10 

15B1-H2-10 

646.80 

15B1-H6-10 

488.70 

Average 

Std.  Dev. 

645.52 

±50.49 

Average 

Std.  Dev. 

447.69 

-34.76 

I5B1-H4-1 

589.60 

15B1-H8-1 

401 .00 

15B1-H4-2 

483.30 

15B1-H8-2 

397.40 

15B1-H4-3 

604.70 

15B1-H8-3 

368.30 

15B1-H4-4 

540.00 

15B1-H8-4 

362.70 

15B1-H4-5 

511  .60 

15B1-H8-5 

428.20 

15B1-H4-6 

539.00 

15B1-H8-6 

382.90 

I5B1-H4-7 

503.90 

15B1-H8-7 

353.20 

15B1-H4-8 

481  .20 

15B1-H8-8 

365.00 

I5B1-H4-9 

539.1  0 

15B1-H8-9 

403.00 

15B1-H4-10 

535.50 

15B1-H8-10 

398.80 

Average 

Std.  Dev. 

532.79 

±40.73 

Average 

Std.  Dev. 

386.05 

±23.50 

k 
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TABLE  XIII 

FAILURE  STRESSES  OF  (:5/0)^  SPECIMENS  WITH  HOLES 


Specimen 

5A1-H2-1 

5A1-H2-2 

5AI-H2-3 

5A1-H2-4 

5A1-H2-5 

Average 

5A1-H4-1 

5A1-H4-2 

5A1-H4-3 

5A1-H4-4 

5A1-H4-5 

Average 

5AI-H6-I 

5AI-H6-2 

5A1-H6-3 

5A1-H6-4 

5A1-H6-5 

Average 

5A1-H8-1 

5A1-H8-2 

5A1-H8-3 

5A1-H8-4 

5AI-H8-5 

Average 


a^(MPa) 


942  ±  41 


836  ±  75 


802  ±  77 


695  ±  34 


V 
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TABLE  XIV 

FAILURE  STRESSES  OF  (=10/0)^  SPECIMENS  WITH  HOLES 


Specimen 

10A1-H2-1 

I0A1-H2-2 

10A1-H2-3 

10A1-H2-4 

I0AI-H2-5 

Average 

10A1-H4-1 

10A1-H4-2 

10A1-H4-3 

10A1-H4-4 

10A1-H4-5 

Average 

10A1-H6-I 
1 0A1-H6-2 
10A1-H6-3 
10A1-H6-4 
10A1-H6-5 

Average 

10A1-H8-1 

I0A1-H8-2 

10A1-H8-3 

10A1-H8-4 

10A1-H8-5 


794 

880 

790 

937 

857 

852  ±  62 

769 

624 

655 

740 

727 

703  r  61 

427 

534 

491 

608 

593 

525  ±  69 

567 
617 
4  91 
458 
478 


Average 


522  1  67 
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TABLE  XV 

FAILURE  STRESSES  OF  (±15/0)^  SPECIMENS  WITH  HOLES 


Specimen 

a^(MPa ) 

15A1-H2-1 

608 

15A1-H2-2 

646 

15A1-H2-3 

800 

15A1-H2-4 

558 

15A1-H2-5 

670 

Average 

656  ±  90 

15A1-H4-1 

412 

15A1-H4-2 

682 

15A1-H4-3 

516 

15A1-H4-4 

600 

15A1-H4-5 

440 

Average 

530  ±  112 

15A1-H6-1 

482 

I5A1-H6-2 

506 

15A1-H6-3 

540 

15A1-H6-4 

527 

15A1-H6-5 

509 

Average 

512  ±  22 

15A1-H8-1 

431 

15A1-H8-2 

363 

15A1-H8-3 

511 

15A1-H8-4 

329 

15A1-H8-5 

432 

Average 

413  ±  70 
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TABLE  XVI 

COMPOSITE 

TOUGHNESS  AND 

m,  SINGULARITY 

LAMINATE 

m** 

[0/=5]^ 

.22 

1300 

[±5/0]^ 

.15 

1140 

[0/±10]^ 

.24 

1060 

[±10/0]^ 

.32 

1260 

[0/±15]^ 

.34 

990 

C±15/0]^ 

.26 

900 

*Dimensions  of  H  are 

MPa 

**FaiTure  stresses  were  modified 
by  a  finite  width  correction. 
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Chapter  I 
INTRODUCTION 

Within  the  past  decade  graphite/epoxy  composite  material 
has  matured  as  a  viable  alternative  to  metals  for  use  in 
high  technology  structures.  Its  comparacively  low  density 
combined  with  high  strength  and  stiffness  make  it  particu¬ 
larly  suited  for  aerospace  applications,  where  weight  reduc¬ 
tion  is  always  a  primary  concern.  Graphite/epoxy's  excel¬ 
lent  performance  as  a  structural  material  stems  from  its 
orthotropic  material  properties.  In  a  strength  designed 
composite  component,  the  fibers  of  individual  plys  can 
ideally  be  oriented  to  produce  a  laminate  with  no  excess 
strength  in  any  direction.  An  equivalent  metallic  component, 
sized  to  carry  a  maximum  load  in  the  direction  of  the  prin¬ 
cipal  stress,  is  usually  overdesigned  when  viewed  from  any 
other  direction.  The  cost  of  this  unavoidable  overdesign,  a 
direct  result  of  metal '  S'-irsotropic  properties,  is  increased 
structural  weight. 

Another  benefit  resulting  from  graphite/epoxy's  ortho- 
tropy  is  the  possibility  'of  creating  specialized  laminate 
properties  by  careful  selection  of  ply  stacking  sequence. 
Aer oe las t i ca 11 y  tailored,  bend ing- twi s t i ng  coupled  forward 
swept  wings  and  near  zero  coefficient  of  thermal  expansion 


space  structures  are  two  of  the  more  sophisticated 
graphite/epoxy  applications  currently  under  study.  The  suc¬ 
cess  of  each  of  these  designs  depends  upon  the  employment  of 
unusual  laminate  properties  which  can  not  be  economically 
duplicated  using  metals  technology. 

Several  aircraft  companies  have  demonstrated  substantial 
manufacturing  cost  reductions  when  replacing  a  metal  compo¬ 
nent  with  a  composite  piece.  These  savings  are  a  result  of 
graphLte/epoxy ' s  molded  fabrication.  Often  up  to  90%  of  a 
billet  or  forging  is  turned  to  scrap  chips  during  the 
machining  of  a  complicated  shape  from  a  metal  block.  In 
comparison,  waste  of  graphite/epoxy  nas  been  shown  to  be  in 
the  10  to  20  percent  range.  Metal  designs  that  for  manufac¬ 
turing  reasons  must  consist  of  several  parts  fastened  by 
rows  of  rivets  can  often  be  produced  as  a  single  composite 
component.  The  reduced  inventory  and  simplified  assembly 
resulting  from  the  use  of  composites  tends  to  lower  manufac¬ 
turing  costs.  As  an  additional  benefit,  the  smaller  number 
of  part  interfaces  and  fasteners  improves  structural  perfor¬ 
mance  by  lowering  the  number  of  stress  concentrations. 

Because  gr aphi te/ epox y  is  a  relatively  new  material, 
insufficient  amounts  of  statistically  significant  data  and 
few  experimentally  proven  analytical  methods  exist.  Without 
these,  the  aggressive  exploitation  of  its  properties  is 
impossible.  Commercial  airframe  manufacturers  have  been 


hesitant  to  make  a  major  commitment  to  the  material  before  a 
more  thorough  understanding  of  its  long  term  behavioral 
characteristics  is  obtained.  In  the  interim  they  have  lim¬ 
ited  its  use  to  easily  removable,  stiffness  designed  compo¬ 
nents  such  as  control  surfaces,  landing  gear  doors  and  fair¬ 
ings.  One  of  the  major  thrusts  in  composite  research,  aimed 
at  providing  information  of  the  type  required  for  structural 
life  predictions,  has  been  the  study  of  graphite/epoxy's 
fatigue  properties. 

Early  fatigue  work  focused  on  determining  the  material's 
response  to  cyclic  tensile  stresses.  Awerbuch  and  Hahn's 
(1)  examination  of  the  behavior  of  unidirectional  coupon 
specimens  is  typical  of  these  tensile  fatigue  investiga¬ 
tions.  In  addition  to  the  usual  measurement  of  life  to  fai¬ 
lure  and  residual  static  strength  over  a  range  of  peak 
stress  levels,  this  study  included  the  effect  of  a  proof 
load  on  the  fatigue  life  distribution.  A  proof  load  is  a 
large  static  load  applied  prior  to  cycling.  It  was  hoped 
that  this  load  would  indicate  some  relationship  between  the 
static  strength  distribution  of  virgin  coupons  and  their 
fatigue  life  by  eliminating  specimens  of  low  static 
strength.  While  their  program  was  not  able  to  conclusively 
show  that  the  strength  anci  life  distributions  of  a  popula¬ 
tion  of  graphite/epoxy  specimens  are  related,  it  did  serve 
as  an  important  source  of  experimental  data  for  statistical 
fatigue  life  prediction  methods. 


Many  of  these  life  prediction  methods  are  based  on  the 
strength-life  equal  rank  assumption  first  introiuced  by  Hahn 
and  Kim  (15).  This  hypothesis  states  that  in  a  group  of 
specimens,  the  one  with  the  greatest  static  strength  will 
also  have  the  longest  fatigue  life.  The  lowest  strength  spe¬ 
cimen  will  be  the  first  to  fail  under  fatigue  loading. 
Static  strength  and  fatigue  life  of  any  particular  specimen 
will  occupy  the  same  position  in  the  two  ordered  sets  of 
strength  and  life  data.  Chou  and  Wang  (5)  compiled  data 
from  various  tensile  fatigue  experiments  and  compared  sev¬ 
eral  strength  degradation  and  failure  models  based  on  the 
equal  rank  assumption.  None  were  found  to  satisfactorily 
describe  every  set  of  experimental  data. 

One  unexpected  result  of  fatigue  tests  using  laminates 
with  circular  holes  and  other  flaws  was  the  discovery  that 
static  strength  often  substantially  increases  during  fatigue 
cycling.  Ramsni  and  Williams  (23),  for  example,  present 
data  displaying  up  to  40%  increases  in  residual  static 
strength  after  5  million  tension-tension  cycles.  Reifsni- 
der ,  Stinchcomb,  and  O'Brien  (31)  explain  this  unusual 
behavior  as  the  simultaneous  interaction  between  a  wear-in 
and  wear-out  process.  Wear-in  occurs  locally  as  the  multi¬ 
ple  damage  modes  available  to  a  composite  laminate  erase  the 
stress  concentration  at  the  discontinuity.  V.' ear-out  con¬ 
sists  of  a  general  degradation  of  the  entire  laminate  by 
repeated  loading.  The  sum  of  these  two  processes  first 
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results  in  a  net  increase  in  strength  due  to  wear-in  and 
later  a  decrease  as  wear-out  becomes  dominant. 

Of  particular  importance  to  the  future  employment  of  gra- 
p'nite/epoxy  are  the  studies  showing  that  specimens  subjected 
to  cyclic  tensile  loads  often  display  fatigue  limit  type 
behavior  when  tested  at  peak  loads  below  approximately  60? 
of  the  laminate's  static  ultimate  strength.  The  fatigue 
limit  is  defined  as  the  load  below  which  a  structure  will 
survive  an  in  finite  (or  at  least  extremely  large)  number  of 
cycles.  Grimes  (11),  in  a  paper  discussing  the  significance 
of  these  findings,  writes,  "This  behavior  is  quite  different 
from  metals,  and  better.  Graphite/epoxy  endurance  limits  are 
not  only  truer  in  a  physical  sense,  but  they  are  higher  in 
magnitude  than  those  artificially  defined  values  for  metals. 
Structures  with  nearly  unlimited  lifetimes  appear  possible, 
at  least  on  paper,  if  the  working  stresses  are  kept  below 
the  initial  facigue  endurance  limit  stress  level." 

The  inclusion  of  compressive  loads  in  the  fatigue  spect¬ 
rum  has  been  shown  to  significantly  reduce  life  to  specimen 
failure.  In  an  often  quoted  report,  Ryder  and  Walker  (37) 
conclusively  demonstrate  this  point  by  presenting  results  of 
tension-tension ,  tension-compression,  and  compr ession-com- 
pression  fatigue  of  two  different  laminates.  The  slope  of 
bot-n  the  tension-compression  and  compr ession- compr ession  S-N 
curves  is  much  steeper  than  the  tension-tension  curve,  with 


the  load  reversal  case  being  the  most  critical.  This  indi¬ 
cates  that  the  introduction  of  compressive  cyclic  loads  can 
greatly  shorten  fatigue  life.  Rosenfeld  and  Kuang  pos¬ 
tulate  that  matrix  strength  is  reduced  to  a  point  where 
fibers  can  buckle  under  a  compressive  load.  This  tends  to 
produce  more  damage  per  cycle  than  an  equivalent  tensile 
load,  and  accounts  for  the  reduction  in  life  compared  to 
tensile  fatigue. 

The  importance  of  compressive  loads  in  graphite/epoxy 
fatigue  is  directly  opposed  to  previous  experience  with 
metals.  This  behavioral  difference  can  be  explained  by  exa¬ 
mining  the  damage  modes  available  to  both.  In  metal  fati¬ 
gue,  crack  growth  perpendicular  to  the  load  direction  is  the 
only  important  form  of  damage.  A  tensile  load,  which  pro¬ 
duces  a  theoretically  infinite  stress  concentration  at  a 
sharp  notch,  causes  extension  by  separating  crack  surfaces. 
Compressive  loads  force  these  opposing  surfaces  together, 
and  therefore  have  little  tendency  to  extend  metallic 
cracks.  Due  to  the  nonhomogeneous  nature  of  a  composite 
laminate,  many  other  damage  modes  besides  crack  extension 
can  exist.  Fiber  breaks,  delamination,  matrix  cracking  and 
fiber/matrix  debonding  are  all  important  damage  modes,  uni¬ 
que  to  composites,  which  may  be  induced  by  static  or  cyclic 
loading.  When  compressive  loads  are  applied,  damage  to  sur¬ 
face  plys  locally  separated  from  the  bulk  of  a  specimen  by 
delamination  can  often  be  accelerated  by  buckling  of  these 
unsupported  regions. 


A  major  source  of  aggravation  to  the  many  composite  lami¬ 
nate  damage  modes  is  the  three-dimensional  nature  of  the 
stress  field  existing  near  free  edges  and  cutouts.  Interla¬ 
minar  shear  along  with  out-of-plane  shear  and  normal  stress 
develop  to  satisfy  the  stress  free  boundary  condition  at 
unloaded  edge  surfaces.  These  stresses,  often  comparable 
in  magnitude  to  the  average  far  field  applied  axial  stress, 
have  been  shown  to  exist  within  a  narrow  edge  region  on  the 
order  of  one  laminate  thick.  Classical  laminate  theory 
ignores  these  out-of-plane  stress  components,  and  therefore 
is  unlikely  to  accurately  predict  the  stress  distribution  at 
precisely  the  location  where  fatigue  is  most  likely  to  begin 
-  at  a  free  edge  or  other  flaw  in  the  laminate.  When  one 
considers  the  inherent  property  scatter,  variety  of  damage 
modes  and  complex  stress  distribution  in  even  a  simple  gra¬ 
phite/epoxy  coupon  specimen,  it  is  understandable  why  past 
analytical,  empirical  and  experimental  attempts  at  fatigue 
life  prediction  have  met  with  limited  success. 

Recent  work  by  Graves  (10)  suggested  a  different  approach 
to  the  investigation  of  composite  fatigue.  Results  of  this 
compressive  fatigue  experiment,  conducted  using  flawed  sand¬ 
wich  beam  specimens  in  four  point  beam  bending,  indicated 

» 

that  pre-failure  damage  propagation  behavior  may  be  more 
consistently  predictable  than  cycles  to  final  failure.  The 
data  showed  that  even  after  develop! ng  extensive  areas  of 
visible  damage,  specimens  could  survive  a  significant  number 


of  additional  fatigue  cycles  before  failure.  While  the 
behavior  of  composite  laminates  after  damage  has  initiated 
is  extremely  important  for  modern  damage  tolerant  aircraft 
structural  designs,  very  little  of  this  type  of  information 
is  available  in  the  current  literature. 

The  research  described  in  this  report  was  conducted  to 
more  rigorously  investigate  the  damage  initiation  and  propa¬ 
gation  behavior  of  graphite/epoxy  laminaoes  with  circular 
holes  ,  during  compressive  fatigue.  Eighty-nine  axially 
loaded,  aluminum  honeycomb  sandwich  specimens  with  gra¬ 
phite/epoxy  face  sheets  of  four  diffent  laminates,  (+45/0)s 
,  (0/+^5)s  ,  (+30/0)s  and  (0/+30)s  were  compressively  cycled 
with  a  minimum  to  maximum  load  ratio  of  R:10.  Fatigue  dam¬ 
age  initiation  and  propagation  were  monitored  on  both  gra¬ 
phite/epoxy  faces  of  each  specimen  using  Moire  interferome¬ 
try,  a  method  of  optically  producing  surface  contour  lines. 
These  patterns  were  photographically  recorded  at  intervals 
for  later  analysis.  Life,  residual  strength,  cycles  to  dam¬ 
age  initiation  and  damage  growth  data  are  reported.  Damage 
initiation  site  and  mode  are  compared  with  the  analytically 
predicted  three  dimensional  stress  field  at  the  edge  of  the 
hole. 
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Chapter  II 
EXPERIMENTAL  METHOD 

2.  1  SPECIMEN  SELECTION 

One  of  the  major  difficulties  of  compressive  testing  is 
the  introduction  of  high  compressive  loads  into  a  specimen 
without  inducing  buckling.  Three  test  method  and  specimen 
configuration  combinations,  designed  to  minimize  buckling 
problems,  were  considered  for  use  in  this  program-  (1) 
externally  supported  coupons,  (2)  sandwich  beam  in  four 
point  bending,  and  (3)  axially  loaded  sandwich  columns. 

A  preliminary  investigation  using  full  and  partial  exter¬ 
nal  coupon  supports  modeled  after  those  described  by  Ryder 
and  Walker  of  Lockheed  (37)  was  conducted.  Resulting  failure 
mode  and  location  seemed  to  be  highly  dependent  on  support 
configuration,  probably  due  tc  the  thinness  (less  than  1mm) 
of  the  six  ply  laminates.  The  further  complication  of  sup¬ 
port  compatibl ity  with  nondestructive  inspection  (HDI)  meth¬ 
ods  led  to  the  rejection  of  supported  coupon  specimens.  Four 
point  bending  sandwich  beams  of  the  type  emoloyed  by  Graves 
(10)  were  also  rejected.  It  was  felt  that  the  concave  curva¬ 
ture  of  the  compressive  face  would  accelerate  the  expected 
deiamination  failure  mode  by  increasing  the  tendency  of  the 
outer  plys  to  buckle  away  from  the  laminate  after  initial 
damage  . 


Axially  loaded  sandwich  columns  were  selected  for  use  in 
this  program  for  several  reasons.  Tne  outer  surfaces  of  the 
specimen  were  free  of  supports  and  thus  accessable  to  NDI , 
while  general  buckling  of  the  2.5nm  thick  sandwich  was  neg¬ 
ligible  at  loads  up  to  the  specimen's  static  failure  load 
(See  Section  4.1).  Both  graphi te/ epox y  faces  simultaneously 
experienced  the  same  load  history;  two  sets  of  initiation 
and  propagation  data  were  tnerefore  generated  with  each  test 
run  . 

2.2  SANDWICH  COLUMN  FABRICATION 

A  total  of  one  hundred  sandwich  beam  specimens,  some  of 
which  were  used  in  the  development  of  test  methods  and  will 
not  be  further  mentioned,  were  fabricated.  Forty-four 
30.5cm.  square,  six-ply  laminates  of  four  different  orienta¬ 
tions  were  constructed  from  30. 5=-.  wide  AS1/5203-5  gra¬ 
phite/epoxy  tape.  Laminates  were  cured  in  batches  of  two  to 
six  using  the  standardized  aut oc 1 ave/ v ac uum  bag/tooling 
procedures  and  cycles  developed  at  MIT  (90).  Cured  laminates 
were  cut  into  five  5x30cm.  coupons  witn  a  water  cooled  dia¬ 
mond  wheel  mounted  on  an  automaci call y  fed  milling  machine. 
The  two  resulting  2.5cm.  wide  edge  strips  were  discarded. 
The  centrally  located  1.25cm.  diameter  hole  in  each  coupon 
was  made  using  a  water  cooled  diam.or.d  drill.  Punch  through 
damage  to  the  bottom  ply  of  one  coupons  during  the  hole 
drilling  operation  was  minimized  by  backing  up  the  piece 
being  drilled  with  graphite/epoxy  scrap. 


The  core  of  the  sandwich  specimens  was  composed  of  three 
separate  pieces  of  6cm.  wide  American  Cyanamid  Aluminum 
Honeycomb,  A  central  17.5cm.  long  72kg/m3  (4.51b/ft3)  test 


section  was  bounded  by  two  6.25cm.  long,  354kg/m3 
(22.11b/ft3)  pieces.  The  heavier  density  core  was  required 
near  the  specimen's  ends  to  withstand  the  2.07  to  3.45  MPa 
(300  to  500  psi)  normal  load  applied  by  the  testing 
machine’s  hydraulic  grips. 

Actual  assembly  of  the  specimens  was  carried  out  in  two 
separate  operations.  Graphite/epoxy  face  sheets  were  bonded 
to  the  three  core  pieces,  which  had  first  been  cleaned  in 
MEK,  using  American  Cyanamid  FM-123  film  adhesive.  Another 
American  Cyanimid  adhesive,  FM-37-2,  was  used  to  bond  the 
individual  core  pieces  during  the  same  operation.  Four 
7.5x6.25cm.  (02/902)5  fiberglass  loading  tabs  were  later 
bonded,  with  a  1.25Gm.  overlap  of  the  core  splice,  to  the 
graphite/epoxy  face  sheets  using  FM-123  (see  Figure  1).  All 
bonding  operations  were  conducted  using  a  hot  press. 

2.3  NONDESTRUCTIVE  INSPECTION  ^  MOIRE  INTERFEROMETRY 

Choice  of  a  suitable  nondestructi ve  inspection  method  was 
of  fundamental  importance  to  the  successful  completion  of 
this  experimental  investigation.  Required  NDI  characteris¬ 
tics  included  in  situ  application,  long  term  stability, 
repeatability,  low  cost,  passive  operation  and  sensitivity 
to  expected  damage  modes.  A  survey  of  methods  indicated  that 
Moire  interferometry  best  met  the  above  criteria. 
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Moire  interferometry  uses  the  interference  pattern  formed 
by  the  relative  motion  of  two  sets  of  fine  parallel  lines  to 
visually  provide  a  map  of  surface  displacement.  Although  not 
actually  used  for  this  experiment,  the  physics  of  this 
method  can  most  easily  be  explained  by  examining  the  more 
commonly  employed  case  of  in-plane  interferometry  for  mea¬ 
surement  of  in-plane  strains.  Typically,  a  grill  of  evenly 
spaced,  parallel  black  lines  and  equal  width  clear  spaces  is 
in  some  manner  fixed  to  the  und e formed  surface  to  be  exa¬ 
mined.-  This  grill  will  defor.Ti  with  the  specimen  surface 
(Figure  2a).  An  identical  reference  grill  applied  to  a 
glass  plate  is  then  placed  near  the  prepared  specimen.  No 
interference  pattern  will  develop  if  the  reference  grill  is 
arranged  to  have  each  of  its  lines  lie  directly  on  top  of  a 
specimen  line  (Figure  2b).  When  the  rigid  reference  glass  is 
moved  the  width  of  one  line  in  a  direction  perpendicular  to 
the  lines,  each  black  line  of  the  reference  grill  will  fall 
directly  above  a  clear  space  on  the  specimen,  and  vice- 
versa.  The  result  is  completely  destructive  interference, 
or  a  uniformly  dark  pattern  (Figure  2c).  By  stretching  the 
specimen,  however,  the  attached  grill  can  be  made  to  contain 
fewer  lines/unit  length  than  the  reference  grill.  This 
creates  a  gradual  change  from  a  totally  dark  to  totally 
clear  patcern.  Each  of  the  resulting  bright  or  dark  lines  is 
termed  a  fringe.  The  distance  bet.-.een  two  equal  intensity 
fringes  is  related  to  the  known  grill  line  spacing;  relative 


displacement  of  points  on  the  surface  can  therefore  be  cal¬ 
culated  by  simply  counting  fringes  between  points  (Figure 
2d)  . 

Delamination  was  the  primary  damage  mode  exhibited  by  the 
composite  specimens  in  this  experiment.  When  the  compressive 
force  applied  during  each  loading  cycle  reached  some  criti¬ 
cal  value,  delaminated  areas  developed  during  fatigue 
cycling  would  buckle.  A  modification  of  the  in-plane  mea¬ 
surement  technique  described  in  the  preceding  paragraph  was 
employed  to  produce  a  visual  indication  of  these  small  out- 
of-plane  surface  deformations  (Figure  3).  To  measure  out- 
of-plane  displacements  with  this  tec hn i que  ,  the  reference 
grill  is  positioned  close  to  and  parallel  with  the  initially 
planar  column  face  to  be  inspected.  A  collimated  light 
source,  placed  at  an  angle  of  between  dO  and  70  degrees  with 
the  surface,  is  used  to  illuminate  the  specimen.  Since  the 
grill  is  located  betv;een  the  light  source  and  specimen,  sha¬ 
dows  of  each  of  its  lines  are  cast  on  the  sandwich  column. 
For  a  reasonably  parallel  beam  of  lignt,  these  shadow  lines 
reproduce  the  reference  grill  pattern  and  take  the  place  of 
the  physically  applied  grill  used  for  in-plane  measurements. 

When  a  portion  of  the  surface  undergoes  out-of-plane  dis¬ 
placement  due  to  local  buckling,  the  shadow  lines  in  that 
region  are  distorted.  This  results  in  the  formation  of 
inoerference  patterns  by  a  process  similar  to  the  one  previ- 


ously  described.  It  can  be  shown  that  the  resulting  pat¬ 
terns  actually  are  surface  contour  lines.  The  difference  in 
height  between  two  successive  fringes  can  be  calculated 
using  : 

AZ  =1/n  tanG 

where  is  change  in  height  between  two  successive  fringes, 
n  is  the  density  of  the  grill  in  lines  per  unit  length  and  0 
is  the  angle  between  the  light  beam  and  the  specimen's  sur¬ 
face  perpendicular. 

In  a  setup  repeated  on  both  sides  of  the  column,  a  TOO 
line/inch  grill  was  mounted  on  a  six  degree  of  freedom 
holder  and  positioned  less  than  0.25cm.  from  the  specimen 
(Figure  4).  Photoelastic,  Inc.'s  Tens-Lac  Type  U-IO-A  Brit¬ 
tle  Lacquer  Undercoating,  a  silver  colored  spray,  was  used 
to  paint  the  sandwich  column  faces  for  increased  Moire  pat¬ 
tern  contrast.  Thirty-five  millimeter  slide  projectors 
angled  50  degrees  from  the  specimen  served  as  light  sources. 
Zoom  lenses  with  knife  edge  slits  mounted  in  front  of  their 
outer  elements  provided  adequately  collimated  beams  (Figure 
5).  From  equation  1,  the  resulting  hz  oetween  adjoining 
fringes  can  be  shown  to  be  0.021cm.  Since  it  was  usually 
possible  to  discern  1/4  fringe  intensity  changes,  out-of¬ 
plane  motion  of  approximately  0.005cm.  could  normally  be 
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One  difficulty  arising  early  in  the  test  program  was  the 
formation  of  an  interference  pattern  between  the  reference 
grill  lines  and  the  peel-ply  textured  surface  of  the  gra- 
phite/epoxy.  It  was  found  that  the  visibility  of  this  sec¬ 
ondary  pattern  could  be  greatly  reduced  by  slightly  rotating 
the  reference  grill  such  that  its  originally  vertical  lines 
w’ere  tilted  approximately  10  degrees.  This  action  produced 
no  visible  change  in  the  observed  damage  patterns. 

2.^  TEST  PROCEDURE 

All  testing  was  conducted  on  a  100,0001b  Material  Test 
Systems  Model  810  testing  machine  operated  in  the  10,0001b 
range.  To  prevent  slippage,  specimens  were  held  in  MTS 
hydraulic  grips  adjusted  to  apply  a  2.07  to  3*^5  MPa  (300  to 
500  psi)  normal  load  to  the  fiberglass  loading  tabs.  The  net 
section  static  compressive  strength  of  all  four  laminate 
types  was  established  using  axially  loaded  columns  identical 
in  design  to  the  fatigue  specimens.  Two  back  to  back  6.35 
mm  3LH  Type  FAE-25-12-S3  strain  gauges,  one  on  each  column 
face  of  most  static  specimens  (Figure  1),  were  monitored 
during  loading  using  BLH  1000  digital  strain  indicators. 
Besides  allowing  the  calculation  of  laminate  modulus  of 
elasticity,  the  linearity  of  the  strain  data  confirmed  that 
no  generalized  specimen  buckling  occurred  before  static  fai¬ 
lure  (Figure  6).  S;.roke  controlled  tests  were  run  at  a 
strain  rate  of  0.0006  per  minute.  When  possible,  load  was 
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released  after  failure  of  just  one  graphite/epoxy  column 
face . 

Fatigue  tests  of  the  axially  loaded  column  specimens  were 
conducted  on  the  same  MTS  testing  machine.  Test  lo^ds  were 
chosen  to  span  the  range  between  static  failure  and  the 
fatigue  limit  with  a  fairly  uniform  distribution  of  data 
points.  Monitoring  of  the  maximum  and  minimum  compressive 
loads  while  cyling  at  20  Hz  and  R=10  revealed  some  diffi¬ 
culty  in  following  the  commanded  sinusoidal  program.  Slight 
reduction  of  the  cyclic  frequency  eliminated  this  system 
resonsance  problem;  thus,  the  majority  of  fatigue  testing 
was  performed  at  19  Hz.  Other  workers  (32)  have  reported 
that  high  cycling  rates  during  a  fatigue  test  can  induce 
temperature  rises  that  alter  matrix  properties.  19  Hz.  was 
chosen  as  a  compromise  between  the  minimization  cf  this 
temperature  effect  and  the  reduction  of  test  time.  In  some 
tests,  frequency  was  gradually  increased  frtm  0.1Hz  to  19Hz 
during  the  first  AOO  cycles  to  reduce  the  likelihood  of 
startup  failures.  -  This  practice  proved  to  dg  unnecessary 
and  was  discontinued  early  in  the  program.  Cycling  was  per¬ 
iodically  halted  for  between  2  and  5  seconds  to  permit  Moire 
pattern  photographs  to  be  taken.  A  static  farce  of  slightly 
more  than  half  the  maximum  compressive  load  was  maintained 
during  this  period.  The  increment  between  successive  pic¬ 
tures  was  varied  to  produce  approximately  e.-en  spacing  on  a 
log(cycles)  plot.  Maximum  exposure  rate  followed  the  sche¬ 
dule  below: 
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CYCLES 

MAXIMUM 

TIME  BETWEEN 

BETWEEN 

NUMBER  OF 

EXPOSURESCHINUTES) 

RA.*;GE(  CYCLES)  EXPOSURES 

EXPOSURES 

AI  li 

0  - 

1000 

200 

6 

0.  175 

1250  - 

4000 

250 

12 

0.219 

A500  - 

10000 

500 

12 

0.439 

11000  - 

20000 

1000 

10 

0.877 

22500  - 

45000 

2500 

10 

2.  193 

50000  - 

120000 

5000 

15 

4  .  336 

1 30000 

-  300000 

10000 

18 

8.772 

320000 

-  600000 

20000 

15 

17.544 

Man  y 

of  the 

low  cycle 

pictures 

were  skipped  during  low 

stress  level  tests  when  experience  indicate!  daoiage  would 
not  initiate  early  in  the  specimen’s  life. 

If  fracture  did  not  occur  by  500,000  cycles,  standard 
procedure  was  to  imi'.ed  ia  tely  conduct  a  residual  strength 
test  by  gradually  increasing  compressive  load  until  static 
failure  was  induced.  Fatigue  loading  of  a  few  specimens  was 
extended,  with  two  specimens  tested  to  one  million  cycles. 
Ten  specimens  were  removed  from  the  testing  machine  before 
final  failure.  Microscopic  examination  of  cross  sections 
taken  near  the  edge  of  the  hole  allowed  comparison  of  Moire 
patcerns  with  actual  laminate  damage. 
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Pseudo- random  specimen  and  test  condition  selection  was 
employed  to  prevent  data  biasing.  Whenever  possible,  speci¬ 
mens  manufactured  from  laminates  cured  in  tne  same  autoclave 
batch  were  distributed  over  a  range  of  load  conditions.  The 
test  program's  history  can  be  reconstructed  from  the  chrono¬ 
logically  sequenced  run  numbers  listed  in  Tables  1  and  £. 

2,5  PHOTOGRAPHIC  TECHNIQUE 

The  intensity  difference  between  light  and  dark  fringes 
of  Moire  patterns  is  normally  low.  Because  a  major  objec¬ 
tive  of  the  investigation  was  the  producion  of  a  photo¬ 
graphic  record  of  damage  initiation  and  propagation,  care 
was  taken  to  increase  the  pattern  contrast  whenever  possi¬ 
ble  .  The  selection  of  a  100  line/inch  grill  was  driven 
mainly  by  this  concern.  Measurement  sensitivity  could  be 
increased  by  employing  a  grill  of  higher  l*ne  density,  but 
only  at  the  expense  of  pattern  visibility. 

Moire  pattern  photographs  were  taken  wish  two  35mm  cam¬ 
eras;  a  motor  driven  Nikon  F2A  equiped  with  a  N'ikon  100mm 
Micro  lens  was  mounted  behind  the  specimen  and  an  Olympus 
OH-1  with  a  Vivitar  90inm  Macro  lens  in  tne  front.  Kodak 
PLUS-X  black  and  white  film  was  used  throughout  the  experi¬ 
ment.  By  intentionally  und  o'"  sxposi  ng  by  2  f/stops  and  dou¬ 
bling  normal  development  time,  contrast  on  the  negative  was 
increased  to  the  point  where  pattern  visioility  was  not  a 
problem  during  measurem.ent  and  printing  of  the  images. 


Chapter  III 
ANALYSIS 

3.  1  THREE  DIMENSIONAL  STRESS  FIELD  PREDICTION 

Experimental  evidence  has  indicated  that  stacking 
sequence  does  affect  both  the  fatigue  life  and  pre-failure 
behavior  of  G/E  laminates.  In  a  more  complete  analysis  of  a 
laminate's  stress  field,  Pipes  and  Pagano  (85)  used  a  finite 
difference  technique  to  show  that  large  interlaminar 
stresses  can  be  developed  in  a  narrow  region  near  the 
straight  free  edge  of  a  coupon  specimen.  Furthermore,  it 
was  found  that  reversing  the  stacking  sequence  of  those  sim¬ 
ple  laminates  changed  the  orientation  of  the  out-of-pi  a/,  c 
normal  stress  from  tension  to  compression.  Since  the  lami¬ 
nate  found  to  have  a  tensile  interlaminar  stress  corresponds 
with  stacking  sequence  experimentally  shown  to  have  a 
shorter  fatigue  life  (22),  it  has  been  hypothesized  that 
tensile  interlaminar  normal  forces  shorten  life  by  acceler¬ 
ating  propagation  of  fatigue  incuced  delaminations. 

A  special  laminated  plate  traction- free  edge  finite  ele¬ 
ment  was  developed  by  Harris  (62)  for  analysis  of  more  com¬ 
plicated  geometries  than  straight  free  edges.  While  this 
element  requires  the  use  of  far  fewer  degrees  of  freedom 
than  would  normally  be  necessry  if  non- spec i al i zed  elements 
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were  employed  to  calculate  the  laminate's  three  dimensional 
stress  field,  the  amount  of  computer  storage  and  time  needed 
to  analyze  even  relatively  thin  laminates  is  still  very 
large.  Also,  the  oscillatory  nature  of  some  of  its  pred¬ 
icted  stress  components  reduces  confidence  in  the  ability  of 
its  forth  order  polynomial  approximation  of  element  stress 
to  accurately  predict  behavior. 

A  more  elegant  method  for  prediction  of  the  three  dimen¬ 
sional  stress  field  in  composite  laminates  is  currently 
under  "development  by  Bar-Yoseph  (52).  By  applying  the  well 
established  fluid  dynamics  technique  of  matched  asymptotic 
expansion  to  the  equilibrium  equations,  the  problem  was 
separated  into  an  inner  region  decribed  by  simple  laminate 
theory  and  an  outer  region  whose  mathematical  description 
included  all  six  stress  components .  Asymptotic  series  des¬ 
criptions  of  tne  stress  field  in  powers  of  tne  small  parame¬ 
ter,  h/1,  the  ply  thickness  nond imen si onal i red  by  some  char¬ 
acteristic  length,  were  substituted  into  the  two  equation 
sets.  The  stress  free  edge  boundary  conditions  were  satis¬ 
fied  exactly  by  imposing  the  condition  that  the  sum  of  the 
inner  and  outer  solutions  of  certain  components  must  be  zero 
at  the  free  edge.  Other  charater istics  of  the  series  were 
determined  from  matching  of  equal  order  terms,  using  the 
fact  that  at  some  intermediate  region  the  solutions  must 
agree.  V/ith  only  terms  of  the  first  order  in  h/1  retained 


for  this  preliminary  application  of  his  analysis,  Bar-Yoseph 
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developed  assumed  expressions  for  the  stress  components  of 
the  boundary  layer  solution  of  symmetric  laminates  in  terms 
of  just  four  exponentially  decaying  position  dependent  con¬ 
stants  per  laminate  ply.  Solution  of  the  3-D  stress  distri- 
outicn  a r our'. 0  a  hole  was  tnerefore  reo^oeo  me  oemer mining 
just  4n-2  constants,  where  n  is  half  the  number  of  plys  in  a 
symmetric  laminate.  Note  that  two  constants  are  determined 
from  symmetry  conditions.  Substitution  of  these  expressions 
into  a  complementary  energy  variational  formulation  resulted 
in  a  set  of  fourth  order  differential  equations.  By  casting 
these  in  the  form  of  an  eigenvalue  problem,  standard  solu¬ 
tion  techniques  could  be  used  to  determine  the  unknown  con¬ 
stants  in  the  stress  functions.  With  these  constants  known, 
the  calculation  of  a  stress  at  any  point  in  the  laminate 
simply  requires  substitution  of  the  coordinates  into  the 
appropriate  equation. 
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Chapter  IV 
DATA  REDUCTIO:; 

A . 1  APPLIED  STRESS  DETERMINATION 

The  axially  loaded  sandwich  column  specimens  used  in  this 
test  program  did  not  buckle  at  loads  below  the  ultimate  com¬ 
pressive  strength  of  the  graphite/epoxy  laminates.  There¬ 
fore,  the  strain  experienced  by  each  of  the  two  face  sheets 
should  be  equal.  If  the  axial  modulus  of  each  laminate  is 
also  assumed  to  be  equal  then  the  applied  stress  can  be 
estimated  by  simply  dividing  the  applied  lead  by  the  total 
of  the  two  face  sheet's  cross  sectional  areas.  With  the 
exception  of  the  far  field  stress  used  in  lam.inate  modulus 
calculations,  all  stresses  in  this  report  are  average  net 
section  stresses  based  on  the  cross  sectional  area  at  the 
hole. 

The  gr a ph i te/ epox y  manuf ac Purer ' s  nominal  thickness  of 
O.SlOlr.m  for  a  6  ply  laminate  was  used  for  all  area  determi¬ 
nations.  Actual  measured  laminate  tnickness  averaged  about 
'iC^  higher.  Microscopic  examination  of  laminates  revealed 
the  source  of  this  discrepancy  to  be  a  thin  textured  layer 
of  epoxy  on  both  outer  laminate  surfaces.  These  layers 
formed  when  excess  resin  material  nardene:  in  the  finely- 
woven  fiberglass  cloth  peel-ply  which  surrounded  the  gra¬ 
phite  /epoxy  during  curing. 
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When  the  peel-ply  was  removed,  an  imprint  of  the  fabric's 
texture  was  left  on  the  laminate's  surface.  The  total  top 
and  bottom  thickness  of  these  resin  layers  almost  exactly 
corresponds  to  the  difference  between  measured  and  nominal 
thickness.  Because  the  resin  rich  surface  layers  are  essen¬ 
tially  unloaded,  it  was  decided  to  ignore  their  contribution 
to  total  laminate  thickness  for  all  stress  calculations  in 
this  report. 

n.2  LAMINATE  MODULUS  DETERHIN'ATION 

Ten  strain  gauged  specimens  were  statically  tested  to 
determine  the  longitudinal  modulus  of  the  four  laminate 
types  used  in  this  program.  Strains,  recorded  every  2001bs, 
during  loading,  were  plotted  versus  the  far  field  stress. 
Modulus  was  calculated  by  finding  the  slope  of  the  best  line 
through  this  data  using  linear  regression  analysis.  The  cor¬ 
relation  coefficient,  an  indication  of  the  linearity  of  the 
30  to  dO  data  points  used  for  modulus  calculaoions,  was  also 
determined  for  each  laminate  (Table  1). 

L.3  DAMAGE  INITIATION  MEASUREMENTS 

Determination  of  the  number  of  cycles  of  fatigue  loading 
endured  by  the  graphite/epoxy  column  faces  before  the  onset 
of  detectable  damage  was  a  primary  objective  of  this  test 
program.  N'egatives  of  the  Moire  patuern  photographs,  taken 
at  intervals  during  fatigue  testing,  were  compared  with  the 


0  cycle  photograph  under  lOx  magnification.  By  starting  from 
a  point  where  damage  was  clearly  visible  and  tracing  back¬ 
wards  through  the  negatives,  it  was  possible  to  locate  the 
earliest  frame  where  damage  could  be  seen.  The  number  of 
cof  respono  ing  uo  Lnis  frame  was  recorceo  as  t,ne 
laminate's  life  to  first  detectable  damage.  Tne  circumfer¬ 
ential  location  of  the  damage  initiation  site  was  found  by 
projection  of  the  frame  containing  first  damage  onto  a  piece 
of  polar  graph  paper.  Positive  angles  were  defined  as  being 
counterclockwise  from  the  load  direction,  with  the  top  of 
the  specimen  as  tested  defined  as  zero. 

Choice  of  the  negative  containing  the  onset  of  fatigue 
damage  was  not  always  obvious,  especially  in  laminates  where 
it  occurred  after  several  hundred  thousand  load  cycles.  In 
some  cases  where  damage  developed  very  slowly,  an  uncer¬ 
tainty  in  cycles  to  initiation  of  one  or  two  negatives  was 
possible.  This  corresponds  to  a  possible  error  of  +■400  to 
40,000  cycles,  depending  on  when  in  the  laminates  life  the 
damage  initiated.  -  This  error  is  relatively  small  when  com¬ 
pared  to  the  order  of  .magnitude  scatter  typical  of  fatigue 
life  d  ata . 

4.4  DAMAGE  AREA  MEASUREHEiJT 

In  order  to  quantify  the  damvCge  propagation  behavior  in 
graphite/epoxy  laminates  after  initiation,  the  surface  area 
bounded  by  the  outermost  Moire  fringe  lines  was  measured  for 
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a  succession  of  individual  frames.  Every  laminate  that 
demonstrated  significant  growth  prior  to  failure  or  test 
termination  was  measured.  Selected  negatives  were  projected 
from  the  rear  onto  a  large  ground  glass  table  top;  the 
resulting  image  was  enlarged  to  approximately  ten  times  life 
size  to  facilitate  measurement.  A  planimeter  was  then  used 
to  mechanically  calculate  the  area  contained  inside  the 
Moire  pattern.  By  also  measuring  the  area  of  the  projected 
hole,  a  simple  ratio  could  be  used  to  determine  the  actual 
damage  area  existing  on  the  specimen. 
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Chapter  V 
RESULTS 


5  .  1  STATIC  TEST  DATA 

Seventeen  specimens,  all  containing  1.25  cm.  diameter 
holes,  were  statically  loaded  to  failure  in  axial  compres¬ 
sion.  Laminate  moduli  and  ultimate  net  section  compressive 
strengths  determinea  from  these  tests  are  listed  in  Tables 
la  and  lb.  Of  particular  note  are  the  high  values  of  R,  the 
linear  correlation  coefficient,  found  for  those  static  tests 
conducted  using  strain  gauges.  Since  R=1  implies  perfectly 
linear  data,  these  numbers  indicate  that  generalized  buck¬ 
ling  of  the  axially  loaded  column  specimens  did  not  occur 
prior  to  compressive  failure  (Figure  5).  Table  1b  shows  an 
unusual  result  of  the  ultimate  strength  measurements .  Alt¬ 
hough  intuition  would  indicate  otherwise,  specimens  with  the 
angle  plys  at  ^5  degrees  were  found  to  fail  ct  higher  aver¬ 
age  compressive  loads  than  those  wit.-,  oh  degree  angle  plys. 
This  observation  has  also  been  rep>:.-ted  by  Graves  (10). 
Since  such  a  limited  amount  of  static  data  is  presented,  it 
is  possible  that  this  effect  is  just  a  manifestation  of  the 
scatter  normally  found  in  composite  properties. 

During  the  tests,  popping  sounds  ceginning  at  loads  lar¬ 


ger 


than  about  QOil  of  the  ultimate  static  load  often  warned 


that  catastrophic  failure  of  the  specimen  was  imminent. 
Fractures  were  usually  perpendicular  to  the  loading  direc¬ 
tion  and  always  through  the  hole.  Because  most  tests  were 
run  in  load  control,  testing  machine  motion  after  failure 
caused  extensive  damage  to  the  fracture  surface.  Testing 
some  specimens  under  stroke  control  eliminated  this  post¬ 
failure  damage  and  resulted  in  some  clean,  sharp  cracks 
across  the  specimen.  This  indicates  delamination  was  prob¬ 
ably  not  a  primary  mode  of  static  failure,  although  Moire 
observation  of  the  surface  of  two  specimens  prior  to  failure 
did  indicate  some  local  delamination  at  the  hole  edge. 

5.2  FATIGUE  LIFE  DATA 

Figure  7,  plotted  from  data  in  Tables  2a  through  2d, 
illustrates  the  fatigue  life  behavior  of  the  four  laminate 
types  investigated  in  this  program.  Results  of  the  static 
tests,  plotted  close  to  the  vertical  axis  and  represented  by 
the  smaller  symbols,  are  also  presented.  A  majority  of 
tests  were  conducted  with  (+m5/0)s  specimens  21  cycled 
to  failure  and  26 'stopped  prior  to  failure.  N'one  of  the 
(4^^5/0)s  specimens  suffered  a  fatigue  failure  when  tested 
at  a  maximum  compressive  load  less  than  300  MPa.  Samples 
tested  at  lower  loads  survived  a  minimum  of  600,000  cycles, 
with  two  tests  extended  to  one  million  cycles  with  no  fai- 
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A  smaller  amount  of  data,  16  points  in  all,  represent  the 
fatigue  failures  of  the  three  other  laminate  types  tested. 
The  data  from  the  (0/+i15)s  fatigue  tests  are  particularly 
interesting.  While  the  static  strength  of  this  layup  was 
found  to  be  nearly  identical  to  the  (+^5/0)s  laminate's,  the 
fatigue  life  was  significantly  lower.  One  specimen  tested  at 
276.1  MPa  failed  after  only  40,000  cycles.  This,  along  with 
the  rest  of  the  data  for  this  laminate,  indicates  that  a 
stacking  sequence  effect  exists  for  the  flawed  laminate  when 
subjected  to  cyclic,  but  not  to  static,  loads. 

No  stacking  sequence  effect  can  be  observed  in  the  lim¬ 
ited  amount  of  (_f30/0)s  and  (0/+30)s  data.  When  tested  at 
the  same  load  levels,  failures  of  both  laminate  types  were 
intermixed  within  a  comparatively  small  life  range. 

5.3  RESIDUAL  STRENGTH  DATA 

The  residual  strengths  for  17  of  the  (+45/0)s  specimens 
which  survived  cyclic  loading  are  listed  in  Table  2b.  These 
data  points,  along  with  the  static  strength  data  for  this 
layup,  are  plotted  in  Figure  8.  Later  sections  of  this 
report  will  show  that  many  of  these  specimens  had  developed 
extensive  areas  of  delaminaticn  prior  to  the  residual 
strength  test.  Even  with  this  visible  damage,  however,  only 
a  slight  reduction  in  load  carrying  ability  from  the  static 


strength  distribution  can  be  seen  after  fatigue. 
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5.^  DAMAGE  INITIATION  DATA 

Measurement  of  life  to  first  observable  fatigue  damage 
was  a  major  objective  of  this  test  program.  The  results  of 
damage  initiation  measurements  are  listed  in  Tables  2a 
through  2d,  with  initiation  data  for  all  four  laminate  types 
consolidated  in  Figure  9.  The  92  points  from  (+45/0)s  tests 
seem  to  demonstrate  a  bilinear  behavior  behavior  on  this 
semi-log  plot.  Points  for  damage  beginning  before  20,000 
cycles  follow  a  negatively  sloping  line  down  to  approxi¬ 
mately  280  MPa,  while  later  initiation  data  is  almost  flat. 
No  damage  was  found  on  any  specimen  tested  below  260  MPa,  in 
two  cases  even  after  1,000,000  cycles. 


Damage  initiation  in  the  (+30/0)s  laminates  falls  within 
the  the  scatter  of  the  (+45/0)s  data.  Both  laminates  with 
the  unidirectional  plys  on  the  outside,  however,  show  signs 
of  damage  significantly  sooner  than  the  average  time  to  ini¬ 
tiation  for  the  cor  respond ing  laminate  with  the  angle  plys 
on  the  outside. 


In  every  specimen  showing  signs  of  damage,  initiation 
occurred  at  the  free  edge  of  the  hole.  Figure  10,  using  data 
from  Tables  2a  through  2d,  plots  the  angular  position  of  the 
initiation  site  measured  counterclockwise  from  the  loading 
direction.  Data  for  each  laminate  type  is  plotted  on  one  of 
the  concentric  circles  representing  a he  hole  edge.  Three  of 
the  laminates,  (0/+30)s  ,(0/+b3)s  ,  and  (+30/0)s  tended  to 


initiate  damage  within  a  10  degree  sector  counterclockwise 


from  the  perpendicular  to  the  load  direction.  The  (+45/0)s 
laminate,  however,  behaved  differently.  In  all  but  a  few 
cases,  initiation  for  this  layup  was  within  +10  degrees  of 
an  intersection  between  lines  +iJ5  degrees  from  the  load 
direction  and  the  hole  edge.  Damage  initiation  behavior  has 
thus  been  experimentally  shown  to  depend  both  on  ply  orien¬ 
tation  and  stacking  sequence. 


5.5  DAMAGE  PROPAGATION  PHOTOGRAPHS 

As  described  in  section  2.3,  Moire  inter ferometry  was 
used  to  observe  damage  initiation  and  propagation  during 
fatigue  cycling.  As  many  as  150  photographs  of  the  Moire 
patterns  on  each  sandwich  column  face  were  taken  as  a  speci¬ 
men  was  tested.  Appendix  A  presents  a  set  of  between  3  and 
16  individual  pictures  selected  to  illustrate  the  progres¬ 
sion  of  damage  for  every  laminate  face  which  showed  a  signi¬ 
ficant  amount  or  growth. 

The  figures  in  the  appendix  have  been  grouped  by  laminate 
type,  and  within  each  laminate  type  are  ordered  chronologi¬ 
cally  as  tested.  Labels  below  each  picture  indicate  the 
number  of  fatigue  cycles  experienced  by  the  specimen  at  the 
time  of  the  photograph.  .The  first  picture  of  a  group  is 
always  at  0  cycles.  It  shows  the  background  pattern  result¬ 
ing  from  a  slight  misalignment  between  the  re  Terence  grill 
and  the  specimen,  as  well  as  some  interference  caused  by  the 


A-34 


interaction  of  the  peel-ply  textured  graphite/epoxy  surface 
and  the  grill  lines.  The  second  picture  is  a  print  from  the 
negative  showing  the  first  sign  of  damage  when  examined  at 
lOx  magnification.  The  last  of  a  series  is  always  the  final 
exposure  made  prior  to  failure  or  the  cessation  of  fatigue 
loading.  In  addition,  some  sets  include  a  few  illustrations 
of  damage  growth  during  statically  applied  loads  of  the 
residual  strength  tests. 

Each  picture  shows  an  area  surrounding  the  hole  in  one  of 
a  specimen’s  graphite/epoxy  face  sheets,  viewed  from  a  point 
along  the  surface  per  pend icular  located  with  its  origin  at 
the  center  of  the  hole.  The  horizontal  dimension  of  the 
print  approximately  covers  the  width  of  the  specimen,  about 
5cm.  The  vertical  dimension,  corresponding  to  the  direction 
of  the  applied  load,  varies  from  set  to  set.  Exact  dimen¬ 
sions  can  be  determined  by  using  the  hole's  known  diameter 
of  1.25cm.  to  scale  any  measurements.  The  test  number  indi¬ 
cated  on  each  figure  corresponds  to  one  same  number  in 
Tables  2a  through  2d,  with  the  letter  designating  which  spe¬ 
cimen  face  is  illustrated.  Figures  with  the  same  number  but 
lettered  a  and  b  are  the  two  sides  of  the  same  column. 

Interpretation  of  the  Moire  patterns  is  simplfied  if  one 
recalls  that  each  fringe  line  shows  an  increment  of  out-of¬ 
plane  motion,  anc  each  picture  represents  a  contour  inap  of 

the  z  direction. 


surface  displacement  in 


Continuous  contour 


lines  indicate  subsurface  delaminations.  Discontinuities  in 
the  contour  lines  occur  at  locations  where  through  cracks 
have  developed.  The  pattern  visible  inside  the  hole  itself 
resulted  from  the  deformation  of  the  underlying  film  adhe¬ 
sive  into  the  core's  honeycomb  cells  during  beam  bonding. 

Inspection  of  the  figures  of  Appendix  A  reveals  several 
interesting  facts.  As  described  in  the  previous  dicussion 
of  Figure  10,  the  damage  initiation  site  is  extremely  con¬ 
sistent  within  a  laminate  type.  Furthermore,  the  pictures 
show  that  each  laminate  type  propagates  damage  in  some  char¬ 
acteristic  pattern  early  in  its  life.  For  the  (+^5/0)s 
laminates,  shown  in  Figures  Ala  through  A41b,  damage  tends 
to  develop  in  a  four  leaf  clover  shape.  Moire  patterns  indi¬ 
cate  that  splitting,  defined  as  through  the  thickness  p'ly 
failures  in  the  fiber  direction,  does  not  usually  occur  for 
some  time  after  initiation.  Contrast  this  behavior  with 
that  of  the  (+30/0)s  laminates,  shown  in  Figures  A44a 
through  A48b.  With  this  layup  the  small  initial  area  of 
delamination,  located  approximately  90  degrees  from  the  load 
direction,  quickly  developed  splits  in  the  direction  of  the 
outer  fibers. 

Tests  of  the  two  laminates  with  0  degree  plys  on  the  out¬ 
side,  Figures  A42a  through  A43b  for  the  (0/-f45)s's  and  Fig¬ 
ures  A49a  through  A54b  for  the  (0/+30)s's  showed  they  ini¬ 
tially  behaved  in  a  manner  similar  to  that  of  the  (+30/0)s 
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layup.  Once  splitting  in  the  load  direction  began,  however, 
the  damage  progressed  much  differently.  Splits  divided  the 
surface  ply  of  each  laminate  into  narrow  regions  which  buck¬ 
led  to  form  small  arches,  each  arch  having  a  length  of 
approximately  one  hole  diameter.  Continued  fatigue  cycling 
caused  the  (0/+30)s  specimens  to  develop  a  stairstep  pattern 
of  these  arches  which  propagated  in  the  direction  of  the 
angle  ply  closest  to  the  surface.  Since  only  a  few  (jfi!5/0)s 
specimens  were  tested,  no  conclusive  statement  regarding 
this  laminates  behavior  after  initiation  can  be  made. 

Graves  (10)  has  also  published  data  describing  damage 
progression  during  compressive  fatigue  of  the  same  'four 
laminates  studied  in  the  current  program.  Surface  damage 
observed  during  periodic  visual  inspections  was  sketched. 
Similarities  between  these  sketches  and  the  Moire  patterns 
of  this  report  can  be  seen,  particularly  during  the  early 
stages  of  development.  Later  in  the  specimen's  life,  how¬ 
ever,  the  results  of  the  two  tests  show  drastically  differ¬ 
ent  behaviors.  Graves  shows  propagation  along  the  specimen 


1 ength  , 

often  ex  ten 

ding  from  one 

end  to  the 

other . 

Damage 

during 

the  current 

investigation 

restricted 

itself 

to  the 

immediate  vicinity  of  the  hole. 

Two  possible  explanations  for  this  difference  are  pro¬ 
posed.  The  sandwich  specimens  used  by  Graves  contained 
0.635cm.  diameter  holes,  half  the  size  of  those  used  in  tnis 
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tified  by  the  sy.r.bol  plotted  at  the  location  of  the  last 
daca  point  measured.  In  most  cases  data  for  both  sides  of  a 
specimen  are  presented.  It  should  be  remembered  that  in 
these  load  controlled  tests,  fracture  of  just  one  of  its 
sides  resulted  in  the  complete  destruction  of  a  specimen. 
When  the  behavior  of  the  two  sides  of  a  specimen  are  com¬ 
pared,  it  often  seems  obvious  which  side  was  responsible  for 
the  final  failure.  The  assumption  that  failure  was  caused 
by  the  side  with  the  most  rapidly  growing  damage  or  the 
greatest  area  was  not  experimentally  tested  for  or  verified. 

Figures  11  through  13  show  data  for  (+45/0)s  specimens 
that  failed  during  facigue  cycling.  Each  figure  contains  the 
data  for  a  small,  arbitrarilly  chosen,  range  of  maximum  com¬ 
pressive  stresses.  Most  of  the  (+45/0)s  curves  can  be 
divided  into  three  regions.  Region  1,  the  nearly  horizontal 
portion,  is  characterized  by  the  initiation  and  slow  growth 
of  damage.  Region  2,  represented  by  a  curved  transitional 
section,  is  an  zone  of  increasingly  rapid  propagation.  The 
end  of  this  region  often  corresponds  to  the  onset  of  major 
splitting  or  the  coalesence  of  two  smaller  damage  areas,  as 
shown  in  Appendix  A.  Region  3,  the  nearly  vertical  section 
of  the  curve,  usually  lasts  for  a  relatively  small  portion 
of  the  specimens  life  but  contains  a  large  percentage  of  the 
total  damage  accumulated  during  fatigue  cycling. 


No  correlation  between  damage  plots  for  the  two  sides  of 
the  same  specimen  can  be  identified.  Sometimes  the  two 
curves  follow  each  other  closely.  At  other  times  damage  on 
one  face  grows  rapidly  while  little  or  none  developes  on  the 
other.  Occasionally  the  two  curves  cross  each  other  one  or 
more  times  before  failure.  The  appearence  of  all  three 
behaviors  suggests  that  no  significant  interaction  between 
the  two  faces  of  a  specimen  occurred,  and  each  pair  of  ini¬ 
tiation  and  propagation  data  can  be  regarded  as  coming  from 
two  separate  tests. 

Plots  which  do  not  follow  the  main  body  of  data,  those  of 
Tests  52b,  58b,  77a,  and  82b  for  example,  begin  Region  3 
type  behavior  early  in  their  lives.  Some  of  these  same 
tests  are  responsible  for  initiation  sites  in  Figure  10  that 
lie  outside  the  four  leaf  clover  pattern.  This  indicates 
that  initiation  in  these  tests  probably  occurred  normally 
but  individual  areao  had  already  joined  by  the  time  a  photo¬ 
graph  was  taken. 

Examination  of  Figures  11  through  13  shows  that  maximum 
damage  area  does  not  always  correlate  with  fatigue  life  or 
maximum  compressive  cyclic  load.  Early  failure  does  not 
necessarily  imply  a  small  final  damage  area  or  a  high  maxi¬ 
mum  load,  as  evidenced  by  test  58b  on  Figure  13*  The  trend, 
though,  is  definitely  toward  larger  amounts  of  damage  devel¬ 
oping  on  the  longest  lasting  specimens,  with  relatively 
small  amounts  of  damage  found  on  those  that  failed  soonest. 
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The  above  descriptions  also  apply  to  the  next  two  graphs, 
Figures  15  and  16,  which  show  damage  progression  in  the 
(+45/0)s  tests  stopped  prior  to  fatigue  failure.  Figure  15 
contains  data  for  specimens  loaded  in  the  maximum  compres¬ 
sive  load  range  which  had  produced  fatigue  failures  in  some 
tests.  Notice  that  several  curves  on  this  plot  indicate 
Region  3  type  behavior  and  suggest  specimens  which  are  prob¬ 
ably  nearing  failure.  In  particular,  note  the  damage  area 
of  Test  63a,  listed  as  10.045cm.  in  Table  2b.  The  pictures 
of  Appendix  A,  Figure  A30a,  show  that  damage  on  this  speci¬ 
men  had  actually  propagated  to  the  left  edge  of  the  coupon 
without  producing  a  fatigue  failure. 

Each  of  the  last  three  plots  of  this  group,  Figure  17 
through  19,  contain  propagation  data  for  both  the  failed  and 
unfailed  specimens  of  one  particular  laminate  type.  Figure 
17,  showing  the  (0/+45)s  data,  indicates  only  one  specimen 
had  developed  sufficient  damage  to  be  plotted.  Others 
either  failed  with  very  small  damage  areas  or  were  removed 
from  testing  soon  after  damage  was  observed.  The  single  plot 
shown  does  indicate  Region  3  behavior  occurred  immediately 
after  initiation. 

The  (+30/0)s  data  in  Figure  18  shows  this  laminate's 
”  » 

behavior  is  very  different  from  the  (+45/0)s's.  After  a 
short  period  of  slow  damage  growth,  rapid  propagation  devel¬ 
ops  with  almost  no  transitional  period.  Fatigue  failures. 
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indicated  by  the  check  marks  next  to  the  individual  run  num¬ 
bers,  occurred  with  relatively  little  observable  damage. 

Contrast  the  (4^20/0)s  behavior  with  that  of  the  (0/+30)s 
specimens  plotted  in  Figure  19-  These  laminates,  with  their 
zero  degree  plys  on  the  surface,  show  the  same  three  regions 
of  damage  propagation  observed  in  the  (+45/0)3  specimens. 
Region  3  growth  is  very  steep,  a  result  of  the  stairstep 
nature  of  the  damage  development  (see  Figures  A49a  through 
A87b) .  Damage  areas  at  failure  were  similar  to  those  of  the 
(+45/0)s's,  and  much  greater  than  damage  found  in  the 
(+30/0)s  tests. 

5.7  MICROSCOPIC  CROSS  SECTION  EXAMINATION 

Specimens  from  ten  fatigue  tests  were  removed  from  the 
testing  machine  with  various  amounts  of  visible  damage,  and 
cross  sections  of  their  laminates  observed  under  250x  magni¬ 
fication.  Three  different  forms  of  damage  were  seen.  Tight 
vertical  cracks  extending  from  one  ply  interface  to  an  adja¬ 
cent  interface,  probably  a  result  of  minor  matrix  failures, 
were  the  most  numerous  form.  Horizontal  spacing  of  these 
vertical  cracks  varied  from  only  a  few  in  the  entire  cross 
section  to  as  many  as  10  per  millimeter  in  areas  of  heavily 
damaged  laminates.  One  laminate  that  showed  no  indication 
of  damage  using  Moire  interferometry  contained  four  of  these 
vertical  cracks  in  a  2cm.  long  cross  section  taken  1mm.  from 


the  hole. 
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A  second  form  of  damage,  generally  corresponding  to  areas 
where  Moire  patterns  were  seen  during  testing,  was  delamina¬ 
tion.  This  type  of  damage  propagated  along  ply  interfaces, 
and  occasionally  jumped  from  one  interface  to  another 
through  a  vertical  crack.  Parallel  growth  of  two  delamina¬ 
tions,  one  above  the  other,  was  only  seen  in  specimens  with 
large  amounts  of  damage.  Major  concentrations  of  the  one 
ply  long  vertical  matrix  cracks  had  at  least  one  of  their 
ends  extending  from  delaminations. 

Through  cracks,  major  intraply  failures  visible  as  splits 
upon  reaching  the  surface,  were  the  least  common  type  of 
damage.  In  all  cases  these  were  associated  with  delaminated 
areas  below  the  surface. 

Damage  in  (+45/0)s  laminates  was  most  concentrated  in  the 
top  two  surface  plys.  Delamination  usually  propagated  along 
ply  interfaces,  but  was  ocassionally  seen  to  drift  into 
intraply  regions  for  short  distances.  In  specimens  with 
large  amounts  of  damage,  delamination  and  vertical  cracks 
did  develop  in  the  plys  near  the  core,  but  damage  to  the 
film  adhesive  layer  between  the  graphite/epoxy  and  honeycomb 
was  never  seen.  It  is  important  to  note  that  some  of  these 
laminates  had  splits  through  the  zero  degree  plys,  indicat- 

I 

ing  that  locally  massive  fiber  ■'eakage  in  the  main  load 
carrying  ply  could  be  sustained  by  a  specimen  without  fai¬ 
lure  . 


Observation  of  laminates  with  their  angle  plys  at  the 
center  of  the  specimen  showed  that  very  shortly  after  initi¬ 
ation,  damage  existed  only  in  these  interior  plys.  A  verti¬ 
cal  split  located  at  the  mininmum  net  section  was  found  in  a 
specimen  that  had  not  yet  developed  the  stairstep  damage  to 
the  surface  zero  degree  plys. 

Figures  20  through  22  contain  50x  micrographs  of  cross 
sections  from  three  different  (+45/0)s  specimens.  The  loca¬ 
tion  of  these  cross  sections  corresponds  to  a  cut  along  a 
vertic’al  line  0.5mm.  to  the  right  of  the  hole.  Because  both 
the  contrast  and  the  magnification  of  these  pictures  are 
low,  much  of  the  minor  damage  observable  during  visual 
inspection  can  not  be  recorded  photographically.  Sketches 
of  observed  damage,  drawn  during  a  slow  ply  by  ply  scan  of 
the  2cm.  long  sections  at  250x,  are  therefore  also  pre¬ 
sented.  Recall  that  while  only  three  samples  will  be  des¬ 
cribed  in  detail,  20  laminates,  each  sectioned  at  two  or 
three  locations,  were  inspected  in  this  manner  to  provide 
the  background  for  the  previous  discussion  of  damage  types. 

Figure  20  shows  damage  sustained  by  the  laminate  of  Fig¬ 
ure  A12a  after  300,000  cycles.  The  numbers  from  -10  to  10 
found  at  the  top  of  the  sketch  indicate  longitudinal  posi¬ 
tion,  measured  in  millimeters,  above  or  below  the  hole  cen¬ 
ter  line.  Many  fine  vertical  matrix  cracks  are  seen  to  have 
developed  in  both  the  inner  and  outer  angle  ply  pairs,  but 


delamination,  indicated  schematically  by  the  heavy  areas 
along  ply  interfaces,  has  only  occurred  in  the  surface  pair. 
The  wide  vertical  lines  in  the  sketch  indicate  splits 
through  the  thickness  of  a  ply,  a  more  severe  form  of  damage 
than  the  much  more  common  matrix  cracks.  Two  such  splits 
were  found  at  this  cross  section.  The  photograph  shows  one 
of  these  splits,  at  location  -5  on  the  sketch,  and  illus¬ 
trates  a  jump  of  delamination  propagation  from  one  interface 
to  another.  A  similar  small  split,  not  yet  through  to  the 
surface,  can  be  seen  forming  at  the  upper  right. 

The  damage  seen  at  this  cross  section  was  relatively 
minor  in  comparison  to  laminates  showing  more  growth  in 
their  Moire  patterns.  The  approximately  1mm.  long  delamina¬ 
tion  located  at  0  in  the  sketch  is  not  visible  in  the  Moire 
picture  of  Figure  A12a,  while  the  2.5mm.  delamination  at  -5 
can  be  associated  with  the  small  spot  on  the  lower  right 
edge  of  the  hole.  This  information  gives  some  indication  of 
the  minimum  delamination  size  detectable  by  the  Moire  inter¬ 
ferometry  method  of  inspection  as  applied  in  this  investiga¬ 
tion. 

Figure  21  is  a  cross  section  of  Test  63b,  Figure  A30b, 
taken  at  1,000,000  cycles.  As  both  the  sketch  and  the  Moire 
pattern  indicate,  much  more  damage  has  developed  close  to 
the  hole  in  this  specimen  than  in  the  previous  example. 
Many  vertical  cracks,  some  as  closely  spaced  as  10  per  mil- 


limeter,  were  found.  Delatnination  between  the  upper  -^<5  and 
0  degree  plys  nearly  spans  the  cross  section.  A  shorter 
delamination  between  the  top  +^5  interface  also  exists.  Ver¬ 
tical  cracks  are  most  numerous  in  the  -U5  degree  ply,  and 
connect  the  parallel  delaminations.  Two  major  splits  have 
propagated  to  the  surface,  but  more  importantly,  a  local 
failure  of  the  zero  degree  plys  has  occurred.  This  split, 
located  at  3  in  the  sketch,  meets  a  delamination  at  the 
inner  -45  ply  interface. 

A  second  cross  section,  taken  approximately  0.5cm.  to  the 
right  of  the  hole  edge,  shows  much  less  damage.  All  forms, 
including  matrix  cracking,  have  been  restricted  to  the  top 
two  plys.  The  two  major  splits  shown  in  the  the  cross  sec¬ 
tion  closest  to  the  hole  are  seen  to  have  propagated  along 
the  fibers  to  this  location. 

The  photograph,  taken  at  the  cross  section  nearest  the 
hole,  shows  the  area  between  10  and  13  in  the  sketch.  A  zero 
degree  fiber  break  is  at  the  far  right,  while  a  major  sur¬ 
face  split  is  at  the  upper  left.  The  complete  delamination 
of  the  top  two  45  degree  plys  is  obvious. 

Figure  22  illustrates  the  laminate  of  Test  44a,  Figure 
A21a.  Moire  patterns  show  damage  in  this  specimen  at 
400,000  cycles  falls  somewhere  between  the  two  previous 
examples.  Microscopic  examination  confirmed  this  observa¬ 
tion.  This  specimen  was  unusual  in  that  a  delamination  runs 


the  length  of  the  cross  section  between  the  inner  0  and  -45 
degree  plys.  More  typical  upper  ply  del aminations ,  along 
with  a  share  of  matrix  cracks  and  splits,  also  exist.  Fai¬ 
lure  in  the  zero  degree  plys  is  again  indicated,  both  in  the 
sketch  and  the  photograph.  Note  that  a  section  of  honeycomb 
core  can  be  seen  in  the  lower  left  of  this  picture. 

To  prevent  any  false  conclusions  from  being  drawn,  it 
should  be  stated  that  of  the  many  cross  sections  examined, 
only  two,  both  taken  very  close  to  the  hole,  contained  zero 
degree  failures.  As  mentioned  above,  difficulty  in  photo¬ 
graphically  recording  the  damage  in  a  cross  section  forced 
pictures  only  of  laminates  with  large,  easily  visible  damage 
areas  to  be  taken.  The  zero  degree  breaks  were  both  inter¬ 
esting  and  obvious,  thus  accounting  for  their  disproportion¬ 
ate  share  of  the  above  damage  photographs. 

5.8  ANALYTICAL  STRESS  FIELD  PREDICTION 

A  modified  version  of  a  computer  program  developed  by 
Bar-Yoseph  (52)  was  used  to  calculate  the  three  dimensional 
stress  field  near  a  circular  hole  in  (+45/0)s  and  (0/+45)s 
laminates.  A  second  program  was  written  to  plot  the  results 
of  the  above  analysis  in  an  easily  in terpretable  form.  It 
had  been  hoped  that  Bar-Yoseph's  analysis,  which  implicitly 
assumes  individual  plys  to  be  composed  of  homogeneous, 
orthotropic  material,  would  be  useful  in  predicting  the 
experimentally  measured  damage  initiation  sites  of  the  lami- 


nates  used  in  this  test  program.  Unfortunately,  as 
observation  of  the  resulting  stress  field  plots  will  show, 
gradients  in  the  region  very  close  to  the  hole  are  extremely 
steep,  with  large  changes  predicted  to  occur  in  a  distance 
corresponding  to  the  diameter  of  a  graphite  fiber  (0.0075 
mm . )  . 

Two  important  points  must  be  recalled  (1)  On  the 
scale  of  this  analysis  graphite/epoxy  is  a  highly  nonhcmoge- 
neous  material,  and  (2)  the  macroscopic  stress  field  near 
the  hole  is  undoubtedly  altered  by  the  microscopic  matrix 
cracks  observed  to  develop  prior  to  and  away  from  the  initi¬ 
ation  sites  of  Figure  10.  Considering  those  facts,  it  is 
not  surprising  that  observed  damage  initiation  sites  did  not 
correspond  to  the  peak  stress  locations  predicted  by  theory. 

The  results  of  this  most  recent  development  in  the  analy¬ 
sis  of  composite  laminate  edge  layers  are  none  the  less  of 
interest  to  researchers  in  several  related  areas,  and  there¬ 
fore  have  been  included  as  Appendix  B.  Figures  B1  through 
B120  are  contour  plots  of  individual  stress  components  cal¬ 
culated  on  various  horizontal  planes.  These  figures  present 
a  top  view  of  each  plane.  Shown  is  l80  degree  sector  of  the 
specimen,  cut  through  the  hole  along  a  line  parallel  to  the 
load  direction.  An  area  from  the  hole  edge  to  a  radial 
location  just  one  ply  thickness  away  has  been  included  in 
the  plot.  Observation  of  Figure  20  shows  that  these  plots 
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span  a  distance  of  less  than  20  fiber  diameters.  The  incre¬ 
ment  between  stress  contours  was  chosen  to  produce  10  cont¬ 
our  levels  per  plot.  Results  are  presented  for  calculations 
performed  at  five  through  the  thickness  locations.  Z=1  and 
1-2  are  ply  interfaces,  while  Z=0.5,  1-5,  and  2.5  are  ply 
mid  planes .  ■ 


The  60  individual  figures  for  each  of  the  two  laminates 
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of  30,  stresses,  nond imen sional i zed  by  the  average  far  field 
axial  laminate  stress,  are  presented  in  the  cylindrical 
coordinates  employed  by  the  analytical  method.  Figure 
titles  indicate  which  stress  component  is  being  plotted, 
with  F  used  in  place  of  the  more  traditional  0  as  an  abbre¬ 
viation  for  stress.  R,  T,  and  Z  of  course  denote  radial. 


tangential  and  out-of-plane  directions.  Note  that  at  the 
hole  edge  the  stress  free  boundary  conditions  have  been 


satisfied  exactly. 


The  second  set  of  30  shows  the  same  data  rotated  into  the 
local  ply  coordinates;  1  is  the  local  fiber  direction,  2  is 
normal  to  the  fibers  and  3  is  out  of  the  laminate  plane. 
FZZ,  which  is  not  changed  by  this  transformation,  is  con¬ 
tained  only  in  the  first  set.  An  additional  plot,  labelled 
as  COMB  in  the  figures,  shows  the  result  of  taking  the 
square  root  of  the  sum  of  the  squares  of  Fi2,  F23,  and  F31 . 


This  value  was  calculated  in  an  unsuccessful  attempt. 


loosely  based  on  the  Tsai-Hill  failure  theory  (67),  to  find 
some  combination  of  stresses  that  would  have  its  maximum  at 
a  damage  initiation  site. 
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Chapter  VI 
CONCLUSIONS 

The  location  of  first  observable  damage  in  four  different 
graphite/epoxy  laminates  has  been  shown  to  be  an  experimen¬ 
tally  repeatable  layup  characteristic.  Most  of  the  92 
(+45/0)s  laminates  that  developed  damage  during  fatigue 
testing  first  showed  signs  of  initiation  in  one  of  four  sym¬ 
metrically  arranged  locations  around  the  hole's  perimeter. 
Damage  in  the  other  three  laminates  began  within  a  10  degree 
sector  measured  counterclockwise  from  the  perpendicular  to 
the  load  direction. 

The  progression  of  the  damage  after  initiation  was  also 
consistent  within  a  laminate  type.  Appendix  A  shows  that 
damage  tended  to  develop  in  patterns  that  can  be  associated 
with  a  particular  layup.  Damage  area  versus  cycle  plots  for 
three  of  the  layups  can  easily  be  distinguished,  one  from 
the  other,  by  their  unique  individual  behaviors  in  the  three 
regions  of  growth  described  by  the  text.  (Insufficient  data 
for  (0/+A5)s  specimens  was  obtained  to  include  them  in  this 
statement)  .  All  three  larpinate  types  showed  a  nearly  hori¬ 
zontal  region  of  damage  growth  (see  Figures  11-19)  early  in 
their  lives.  The  (  +  ^15/0)5  specimens  tended  to  gradually 
increase  their  rate  of  damage  growth  in  a  second  region. 
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Tnis  transition  zone  was  much  less  visible  on  the  plots  of 
the  two  laminate  types  with  30  degree  angle  plys.  A  third 
period  in  the  damage  progression  for  all  three  laminates  was 
the  region  of  very  rapid  propagation  snortly  before  failure. 
The  damage  area  at  failure  of  the  (+30/0)s's  was  signifi¬ 
cantly  smaller  than  that  of  the  other  two  laminate  types. 

An  important  experimental  result  was  the  lack  of  correla¬ 
tion,  in  all  but  the  broadest  sense,  between  final  damage 
area  and  either  maximum  fatigue  load  or  life.  '  In  all  four 
laminates  tested,  a  majority  of  the  compressive  load  was 
carried  by  the  zero  degree  plys.  Post-fatigue  inspection 
showed  damage  concentrated  mainly  in  the  angle  plys.  This 
suggests  that  better  correlation  between  damage  area  and 
load  may  be  found  in  tests  of  laminates  that  initiate  and 
propagate  damage  in  their  load  carrying  plys. 

Many  specimens  showed  a  rapid  growth  of  damage  area 
immediately  prior  to  failure.  This  result  suggests  that  some 
warning  of  imminent  specimen  fracture  may  be  obtained  by  the 
periodic  monitoring  of  fatigue  damage,  with  emphasis  placed 
on  the  rate  of  change  of  damage  growth  rather  than  its  phy¬ 
sical  si  ze  , 

.  forty  (+45/0)s  cross  sections,  taken  at  different 
•  "  r.::c  respect  to  the  hole  and  from  laminates  show- 
•  .  1  y  of  damage  in  their  Moire  patterns,  have 

A  ;  tie  sequence  of  damage  propagation 


for  this  type  of  laminate  has  been  developed  from  these 
observations.  Initial  damage  in  a  ply  takes  the  form  of 
vertical  matrix  cracks  which  originate  either  from  high 
static  loads  (first  cycle  damage)  or  fatigue  cycling.  These 
cracks,  which  almost  always  stop  at  ply  interfaces,  act  as 
nucleation  sites  for  delamination  by  perturbing  the  already 
large  interlaminar  stresses  that  occur  at  these  locations. 
Delaminations  then  grow  to  a  size  where  a  ply,  weakened  by 
the  combination  of  wear-out  type  fatigue  damage  and  lack  of 
support  from  adjoining  plys,  can  no  longer  carry  the  maximum 
applied  load  and  splits  along  the  fibers.  The  load  carried 
in  the  split  ply  must  be  transfered  to  other  lamina,  further 
accelerating  damage  propagation  and  producing  the  Region  3 
type  behavior  described  in  the  previous  section. 

Moire  interferometry  proved  to  be  a  simple,  quick  and 
effective  method  for  the  detection  of  subsurface  damage  in  a 
composite  material.  Microscopic  examination  of  damaged  spe¬ 
cimens  showed  that  deiaminations  with  a  length  of  approxi¬ 
mately  2mm.  could  be  routinely  detected  in  the  top  few  plys 
of  the  laminate.  The  method  was  less  useful  in  the  detec¬ 
tion  of  cracks  and  splits,  since  these  forms  of  damage  do 
not  tend  to  produce  out-of-plane  thickness  changes.  Since 
the  patterns  produced  are- visible  to  the  naked  eye,  real 
time  observation  of  damage  propagation  is  possible.  This 
characteristic  of  Moire  interferometry  gives  it  an  obvious 
advantage  over  many  of  the  more  commonly  employed  composite 
inspection  methods,  especially  in  experiments  of  this  type. 


An  analysis  of  tu-i  three  dimensional  stress  field  near 
the  hole  was  unable  to  predict  the  damage  initiation  sites 
of  two  of  the  laminates  tested.  Because  graphite/epoxy  is 
highly  nonhomog en eous  on  the  scale  of  the  analysis,  its  fai¬ 
lure  was  not  surprising.  It  seems  that  an  even  more  compli¬ 
cated  approach,  taking  into  account  the  material's  inherent 
inhomogeneity,  will  be  required  to  accurately  predict  initi¬ 
ation  sites  in  an  arbitrarily  stacked  laminate. 
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